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Abstract 

We perform a fully relativistic calculation of the 2 H(e, e'p)n reaction in the impulse approxima- 
tion employing the Gross equation to describe the deuteron ground state, and we use the SAID 
parametrization of the full NN scattering amplitude to describe the final state interactions (FSIs). 
The formalism for treating target polarization with arbitrary polarization axes is discussed, and 
general properties of some asymmetries are derived from it. We show results for momentum distri- 
butions and angular distributions of various asymmetries that can only be accessed with polarized 
targets. 
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I. INTRODUCTION 



There are many interesting questions to be answered in investigating exclusive electron 
scattering from the deuteron: what does the nuclear ground state look like at short distances, 
are there any six-quark contributions to the wave function, when does a description in terms 
of hadronic degrees of freedom break down? In order to answer any of these questions, 
a precise understanding of the reaction mechanism is mandatory. Final state interactions 
are the most relevant component of the reaction mechanism at GeV energies, but meson 
exchange currents and isobar states will also contribute. The fact that the deuteron is the 
simplest nucleus enables us to study all facets of the reaction mechanism in great detail. 
Anything that can be gleaned from the deuteron will be highly useful for heavier nuclei. 
Exclusive electron scattering from nuclei is one type of reaction where one may observe color 
transparency and the deuteron itself provides a laboratory for the study of neutrons, 
e.g. the neutron magnetic form factor [2]. The short range structures studied in exclusive 
electron scattering might even reveal information about the properties of neutron stars 
For some recent reviews of exclusive electron scattering, see e.g. 

Recently |z|, we performed a fully relativistic calculation of the D(e, e'p)n reaction, using 
a relativistic wave function js| and NN scattering data ^ for our calculation of the full, 
spin-dependent final state interactions (FSIs). The main difference to many other high qual- 
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12] or a diagrammatic 



ity captions ns,ng the generalised eikona. approximation 
approach [13| is the inclusion of all the spin-dependent pieces in the nucleon-nucleon ampli- 
tude. Full FSIs have recently been included in fl4|. Several experiments with unpolarized 
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deuterons are currently under analysis or have been published recently, 
There are also new proposals for D(e,e'p) experiments at Jefferson Lab 19]. 

In we focused on observables that are accessible for an unpolarized target and an 
unpolarized nucleon detected in the final state. The spin-dependent pieces in our FSI cal- 
culation were particularly relevant for the fifth response function, an observable that can be 
measured only with polarized electron beams. Naturally, experiments with polarization of 
the target or ejectile are harder to perform than their unpolarized counterparts. However, the 
extra effort allows one to study otherwise inaccessible observables that are rather sensitive 
to certain properties of the nuclear ground state and the reaction mechanism. In this pa- 



per, we investigate the asymmetries that can be measured with a polarized deuteron target. 
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These observables are of particular interest to us as we have a precise, fully spin-dependent 
description of the final state interactions. As before, the focus of our numerical calculations 
is the kinematic region accessible at GeV energies, i.e. the kinematic range of Jefferson Lab. 
Currently, some deuteron target polarization data that were taken in Jefferson Lab's Hall 
B are being analyzed 20[]. At lower energies, measurements of the beam- vector asymmetry 
A^ d have been performed at NIKHEF 21] and at MIT Bates , |22). A formalism was 
developed within a non-relativistic framework, and calculations of various asymmetries at 
lower energies were performed in 23j. The tensor asymmetry has been discussed at higher 
energies within a Glauber theory approach, with just a central FSI, in [24]. 

The paper is organized as follows: first, we establish the general formalism necessary 
to calculate response functions for polarized targets, and we discuss how to perform these 
calculations in two different coordinate systems. Then, we continue with the calculation of 
asymmetries, and with the issues presented in using the experimental convention of mea- 
suring polarizations along the electron beam direction versus the theoretical choice of po- 
larization axis along the three-momentum transfer q. In the next section, we present our 
numerical results, in a kinematic region relevant to experiments at Jefferson Lab. We show 
both momentum distributions and angular distributions, and we discuss the contributions 
of the various spin-dependent parts of the final state interactions, as well as the influence of 
the ground state wave function. We conclude with a brief summary. 



II. FORMALISM 



A. Differential Cross Section 



The standard coordinate systems used to describe the D(e,e'p) reaction are shown in 
Fig{TJ The initial and final electron momenta k and k' define the electron scattering plane 
and the xyz-coordinate system is defined such that the z axis, the quantization axis, lies 
along the momentum of the virtual photon q with the x-axis in the electron scattering plane 
and the y-axis perpendicular to the plane. The momentum p of the outgoing proton is in 
general not in this plane and is located relative to the xyz system by the polar angle 6 P and 
the azimuthal angle <j) p . A second coordinate system x'y'z' is chosen such that the z'-axis is 
parallel to the z-axis and the x'-axis lies in the plane formed by p and q and the y'-axis is 
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FIG. 1: (Color online) Coordinate systems for the D(e,e'p) reaction, k and k' are the initial 
and final electron four-momenta, q is the four-momentum of the virtual photon and p is the four- 
momentum of the final-state proton. 



normal to this plane. 



The general form of the D(e, e'p) cross section can be written in the lab frame as 26|, |27] 



da 5 



del 'dVt e dVt p 7 , 



m p m n p p 



&Mott free 



87T 3 M d 

~\~h vlt'Rlt 1 + hvx'Rr 



v L R L + v t Rt + vttRtt + v lt Rlt 



(1) 



where Mj, m p and m n are the masses of the deuteron, proton and neutron, p p = p± and Q p 
are the momentum and solid angle of the ejected proton, e' is the energy of the detected 
electron and Q e is its solid angle, with h — ±1 for positive and negative electron helicity. 
The Mott cross section is 



O'Mott 



a cos(6' e /2) 
2e sin 2 (# e /2) 



and the recoil factor is given by 

fr 



up p - E p q cos 6 P 
M d p p 



(2) 



(3) 
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The leptonic coefficients vk are 



v T = + tan — (5) 



2q 2 ' ' " 2 

Q 2 

2q 



V TT = ( 6 ) 



^ = - W ^ + tan- (7) 

Q 2 $e / n \ 

Vlt' = f= — tan — (8) 

Qe / Q 2 o /r>\ 

ur/ = tanyW — + tan z — (9) 

Within this general framework, we have two options for evaluating the response functions: 
first, we will give expressions for the response functions in terms of matrix elements that are 
defined with respect to the electron plane, i.e. the xyz plane. These matrix elements are 
implicitly dependent on <f> p , the angle between hadron plane and electron plane, and these 



are the responses used e.g. in [26|]. Second, we give expressions for the responses in the 
x'y'z' plane. All quantities given relative to the x'y'z' coordinate system are denoted by a 
line over the quantity. The current matrix elements, and therefore the response functions, 
in the x'y'z' coordinate system do not have any <p p dependence. It is much more practical 
to evaluate the responses in the x'y'z' coordinate system. The commonly used responses in 
the xyz system can then easily be obtained by accounting for the (j) p dependence explicitly, 
see eq. (fT7|) below, instead of newly evaluating matrix elements for each value of <J) p . 

Note that both coordinate systems use the same quantization axis: the z axis and the z' 
axis are parallel. We will discuss using a different polarization along the beam, as commonly 
done by experimentalists, in the next subsection on asymmetries. 

The hadronic tensor for scattering from polarized deuterons is defined as 



W XirXj (D)= (Pl S l,P2S2-,(-)\J\i / \PX d )*(PlSl,P2S2](-)\J\ J \P^d}p\ d \> d (10) 

where 

J ±1 = T^=(J 1 ± J 2 ) (11) 
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and 

Jo = J° (12) 

is the charge operator. The notation (— ) in the final state indicates that the state satisfies 
the boundary conditions appropriate for an "out" state. The deuteron density matrix in the 
xyz-fr&me is 

( 1 + yfl T 10 + ^T 20 - y^TA + T* x ) VST* 2 ^ 

-y/i(T n + T 21 ) l-V2T 20 -y/ii^-T^) (13) 
V -\[l{Tii - T 21 ) 1 - y/fTxo + ^T 20 

and the set of tensor polarization coefficients is defined as 
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D = {U,T 10 ,T 11 ,T 2O ,T 21 ,T 2Z } (14) 

with U designating the contribution from unpolarized deuterons. The derivation of the 
density matrix and the conventions used are described in the Appendix. 
The response functions in the xyz-frame are given by 

R L (D) = w 00 (D) 

R T {D) = w u (D) + w-i-i(£>) 
R TT (D) = 2^(w 1 _ 1 (D)) 
R LT (D) = -2% 01 (D)-«; _ 1 (D)) 
Rlt'(D) = -2^(w i(D) + wo-i(D)) 

R T ,{D) = w 11 {D)-w. 1 . 1 {D) (15) 

Now we proceed to write down expressions for the responses in the x'y'z' coordinate 
system. Calculating the responses in this system offers a faster alternative to the above 
calculation, which requires a new evaluation of the current matrix elements for each <p p 
value. The response functions defined above are implicitly dependent upon the angle (f> p 
between the electron plane and the hadron plane containing the proton and neutron in the 
final state. This dependence can be made explicit by noting that 

(pisi;P2s 2 ; (-)l Jx,, \P\ d ) = e *(A d +A 7 - sl - S2 )^ (piSi;p2S2 . ( _ } | lPXd) (16) 

where the line over the matrix elements is used to indicate that they are quantized relative 
to the x'y'z' coordinate system. The hadronic tensor can then be written as 

w x ,, Xj (D) = e-^'- x ^w x ,, Xl (D) (17) 
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where 

w x ^(D) = Yl (Pi s uP2S2](-)\Jx'jPy d )*(piSi,P2S2]( y -)\Jx^\P\d)Px d \' d (18) 

and 

ft* = e^-^p D XdX , d (19) 

is the density matrix defined relative to the x'y'z' coordinate system. 

Using eq. (TTTT) and the definition of the responses in the xyz system, eq. (EES]), the 
response functions in the x'y'z' system then become 



Rl(D) 


= R [ 1\d) 




R T (D) 


= R { t ] (D) 




Rtt(D) 


= Rtt(D) cos 20 p + R^ (D) sin 




Rlt(D) 


= R { l T (D) cos (p p + R { It (D) sin q 


'p 


Rlt>(D) 


= r[ t , (D) sin (f) p + R^It' {D) cos < 




R t >(D) 


= R^P(D) 





where the reduced response functions for the two classes I and II are defined in terms of the 
hadronic tensors as 







^sf = 


Wqo(D) 








tZJi j i(D)+iF_i > _i(D) 


Rjiji(D) 


i 

= -R rT ( 


^)t1 7) = 


2»(TU a ,_ 1 CD)) 


Rj"ip (D) 


i 

= Rrprp 




= 23(u; 1 ,_ 1 (D)) 


Rl ] t(D) 


i 




-2&(Woi(5)-Wo_i(5)) 


Rl?(D) 


i 

= y^-^LT 

i 




= 23f(u; 01 (D)+TZ;o_i(S)) 


RlUd) 


i 


irPyrr = 


= 29f(Woi(S)-W - 1 (S)) 


R { It\d) 


= Rlt'i 




= -2K(W 01 (D)+W _ 1 (D)) 




i 

= y^-R T / 
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= Wi,i(S)- : uJ_i I _ 1 (D) J 



where 



T\ ' G {U,^(T U ),T 20 ^(T 21 )^(T 22 )} 
t[ H) G {T 10 ,$t(T n ),%(T 21 ),Z(T 22 )} (22) 



and 

.(-0 r- f-t „ rJl jst' 



T; 



T; 



{II) 



^ t 1 ; T ll> T 20, T 21 , T 22 ) 

Gjno,^,^,^}. (23) 



The r-matrices are defined by flA20j) . (IA28I) and flA29j) . The type / and // response functions 



can be obtained directly by noting that the density matrix can be written as 

?=5( 1 + E^." , + E-. ( " , r.'" ; ] ■ (24) 

\ i i 

Defining a set of projected hadronic tensors as 



^,A 7 (rf' /7) ) = \ (P^P^; (-)l A> |PA^)*(p lSl ;p 2 5 2 ; R| J Xi \PX d )(r^ n) 
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(25) 

the type / and II response functions are then obtained by replacing the hadronic tensors 
on the right-hand side of the expressions in (j2ip with each appropriate projected hadronic 
tensor in turn. Note that the r matrices satisfy 



and 



t >" } ) > „ = L • < 27 > 

/ — A — A V / AA 



B. Symmetries of the Current Matrix Elements 

The current matrix elements used here are defined in j3] . The matrix elements quantized 
in the hadron plane x'y'z' can be shown to satisfy the symmetry 



(p lSl ;p 2 s 2 ; (-)| J A7 \PX d ) = (-l^+V-i-^ - Sl] p 2 - S2] (_)| J A _ 7 |P - A d ) (28) 
by starting with 

iZ 2l u(p,s) = (-l) 1 * +s u(p,-s) (29) 



which relies on the fact that the nucleoli momenta have, by construction, no y' component 
when quantized in the hadron plane. 

Application of parity and time reversal to these matrix elements requires that 

(pi8i;P28 2 \ HI Jx, \PXd) = (-l) x ^- s ^ (P - X d \ J A _ 7 \ Pl - Sl ;p 2 - s 2 ; (+)) . (30) 

Combining this with ( |28i) gives 

{pis 1 ;p 2 s 2 ;(-)\J\ 1 \P\d} = (PXd\J\ 1 \piSx;p 2 s 2 ;(+)) 

= (PiSi,p 2 s 2 ;(+)\JxjPX d )* . (31) 

In the plane-wave approximation there is no difference between the (— ) and (+) boundary 
conditions. So in this approximation the current matrix elements are real. 

C. Asymmetries 

The simple form of ([T]) is due to the choice of quantization axis associated with the plane 
determined by the virtual photon momentum and the ejectile momentum. In practice, 
the polarization coefficients are determined relative to a coordinate system fixed in the 
laboratory with the axis of quantization along the electron beam momentum. This can be 
easily accommodated by rotating the density matrix. The relationship between the density 
matrix in the x'y'z' coordinate system and the system with the quantization axis z" along 
the electron momentum k and with y" parallel to y is 

P^' d = E D Ia(-^ ^ °) D k A'(-&» °)Paa' ( 32 ) 

AA' 

where the tilde denotes the density matrix for the x"y" z" coordinate system and 9k q is the 
angle between the beam momentum k and the momentum transfer q. The polarization 
coefficients Tjm can be found as functions of the Tjm by using 

T J0 = Tr(Tj p) 
U(T JM ) = ^Tr 



(33) 

The response functions for the x"y"z" coordinate system can by found by using these in 
«2U). 
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SR t— 
r JM P 



Q(T JM ) = -Tr 



9 t- 
TjM P 



The asymmetries that we will calculate here involve the case where T 10 is nonzero with 
all other polarization coefficients equal to zero, or where T 20 is nonzero with all other polar- 
ization coefficients equal to zero. In the first case, 

T w = cos 6 kq f 10 
3?(T U ) = sin 6> fcg cos p Tio 

3(Tn) = -j= sin 9 kq sin 4> p f w 

T 2M = (34) 

while in the second case 

T 1M = 

T 20 = -(l + 3cos20 fcg )T 2O 

— I~3 

$(T 2 i) = ~\ - sin 26> fcg cos P T 2O 
V o 

— I~3 

S(T 2 i) = \ - sin 26 kq sin <f) p T 20 
V o 

3?(T 22 ) = y^(l-cos2^)cos20 p T 2O 

S(T 22 ) = --y^(l - cos26» fcg ) sin20pT 2O . (35) 
A similar relation between the xyz and x"y"z" coordinates systems is given by 

P^y d = ^ d l d A(0kg)d\, jA/ (9 kq )p® A , . (36) 

AA' 

Then, 

Tjm = Tjm(D) . (37) 

The relations between the polarization coefficients can be obtained from (1341) and (I3l)|) by 
setting cf) p = and making the replacements T JM — > T JM . 
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The single and double asymmetries for these two polarizations are defined as 







v L R L (T 10 ) + v T R T {T w ) + v TT R TT (T w ) + v L tRlt(T 10 ) 

T 10 S 

v L R L (T 20 ) + v T R T (T 2 o) + v T tRtt(T2q) + vltRlt(T 20 ) 

T 2 o S 

vlt>Rlt'(Tiq) + Vt'Rt'(Tiq) 
vlt'Rlt'(T 2 q) + vt'Rt'(T 20 ) 



(38) 



T 2 o^ 



where 



£ = v L R L (U) + v T R T {U) + v TT R TT {U) + v LT R LT (U) . 



(39) 



Here Ri(T 1Q ) and Ri(T 20 ) denote the response functions where only T 10 is nonzero or only 



T20 is nonzero. Ri(U) denotes the unpolarized response functions. 

Using the definitions of the asymmetries, the expressions for the Tjm as a function of 
the T and the definitions of the response functions in the x'y'z' system, one obtains the 
following symmetry relations with respect to 4> p : 



III. RESULTS 

All results are shown for a quantization axis along the beam direction, as measured in 
experiments, not along the direction of the three-momentum transfer. 

1. Momentum Distributions 

In Fig. [21 we show the four asymmetries for a four-momentum transfer of Q 2 = 2 GeV 2 
and x = 1. These kinematics correspond to quasi-elastic scattering. Note that in the 
plane-wave approximation the asymmetries and A T ed vanish. They are non-zero only 




-^(360° 
Aj(360° - <j> p ) 
i4£(360° - cj> p ) 




(40) 
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FIG. 2: (Color online) The asymmetries A d (panel (a)), A T d (panel (b)), A\ d (panel (c)) and A^ d 
(panel (d)) for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , x = 1 GeV, and (j) p = 35° are shown 
calculated in PWIA (dotted line), with on-shell FSI (solid line), and with on-shell and off-shell FSI 
(dashed line), as a function of the missing momentum. 

when the FSIs are included. This is predicted in non-relativistic PWIA calculations, and 
our relativistic approach does not change this feature. In the cases where the asymmetries 
are non-zero for PWIA, the inclusion of FSIs leads to a shift, and a slight distortion, of 
the features that are already present in the asymmetries. The dips and bumps become 
narrower when FSIs are included, and they appear at somewhat lower missing momenta. The 
difference between just on-shell FSI and full FSIs including on-shell and off-shell distributions 
is very small. The largest off-shell FSI effects are present for larger missing momenta in A T ed . 



From our discussion in section III B[ we can now explain the observed behavior of the 
asymmetries in PWIA: all PWIA current matrix elements are real, and so any response 
that consists of taking the imaginary part of any part of the hadronic tensor will vanish 
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in PWIA. When consulting eq. (1221) . we see that the vector asymmetries with Tio 7^ are 
associated with the class II responses, and that the tensor asymmetries with T20 7^ are 
associated with the class I responses. From its definition, we can see that A v d is associated 
with the vector, i.e. class II, contributions to the L, T, TT, and LT responses. The L and 
T responses have no class II versions, and the class II versions of the TT and LT responses 
are proportional to the imaginary part of certain pieces of the hadronic tensor. Thus, A% 
vanishes in PWIA. A similar argument shows that Aj d must vanish in PWIA, whereas the 
other two asymmetries will always have non-zero contributions. This argument was made in 
the hadron plane, in the x'y'z' frame. It is also valid when the quantization axis is rotated, 
as the rotation itself will not lead to a non-zero value for an asymmetry that vanishes for 
one set of quantization axes. 
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FIG. 3: (Color online) The asymmetries A d (panel (a)), A^ (panel (b)), A^ d (panel (c)) and A^ d 
(panel (d)) for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , x = 1.3 GeV, and <p p = 35° are shown 
calculated in PWIA (dotted line), with on-shell FSI (solid line), and with on-shell and off-shell FSI 
(dashed line), as a function of the missing momentum. 
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In Fig. [3j we show the four asymmetries for a four-momentum transfer of Q 2 = 2 GeV 2 
and x = 1.3. These kinematics are away from the quasi-elastic peak, and we expect off- 
shell contributions to the FSIs to be more relevant here. We have observed the increase 
in relative importance of the off-shell FSIs already for unpolarized observables in Due 
to the chosen kinematics, smaller values of the missing momentum are not accessible. As 
for the quasi-elastic kinematics shown above, A\ and A T edl are non-zero only once FSIs are 
included, and the FSIs shift the bumps and dips to lower missing momenta. The shift to 
lower momenta is much smaller here than for the quasi-elastic case, though. In contrast to 
the x = 1 kinematics, the off-shell FSIs now play a more prominent role. The differences 
between just on-shell FSIs and off-shell and on-shell FSIs are large for A d and Aj d , and 
they are apparent already at low missing momentum. For the two other asymmetries, 
A^ d and A J, the differences are less pronounced and are most significant at the largest 
missing momenta considered here. Having a non-zero asymmetry already in PWIA makes 
the asymmetry less sensitive to off-shell effects: if the PWIA results are non-zero, the FSIs 
are very relevant corrections, and the off-shell FSIs are less significant corrections of the 
correction; if the PWIA results are zero, the FSIs provide the entire asymmetry, and the 
off-shell FSI corrections are relevant. 

The asymmetries that we investigate here also have a dependence on the azimuthal angle 
cf) p of the outgoing proton. The two sets of figures above were shown for a value of <p p = 35°. 
This value was chosen to avoid any special cases for p = 0°, 45°, or 90°. However, the 
overall 4> p dependence is interesting, too. We show this dependence for all four asymmetries 
in Fig. H]in a three-dimensional plot. 

One can see that for A d , the broad bump and dip structures observed for (fi p = 35° turn 
into a broad dip and bump for <ft p values above 180°, inverting the original, low <p p structure. 
A very similar inversion of the structures is observed for A T ed : the broad ridge at lower (f) p 
turns into a valley for large (fi p , and the sharp dip at low missing momenta and medium 
<p p turns into a peak at <p p > 180° For the other two asymmetries, A T d and A^ d , the plots 
are symmetric around <ft p = 180°. This is the behavior predicted by eq. (HOT) for the four 
asymmetries. 

For the kinematics away from the quasi-elastic peak, for x = 1.3, the same type of 4> p 
dependence and the same <p p symmetries are observed, and we therefore do not display a 
separate figure. The asymmetries reach much larger maximum values for x = 1.3, though. 
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FIG. 4: (Color online) The asymmetries A d (a), A d (b), A^ d (c), and A T ed (d), for a beam energy 
of 5.5 GeV, Q 2 = 2 GeV 2 , and x = 1 are shown calculated with on-shell FSI as a function of the 
missing momentum and the proton's azimuthal angle. 

2. Angular Distributions 

We now discuss our results for angular distributions. Note that for the FSI calculations, 
there is a limit to the kinematic region we can calculate for, as the proton-neutron scattering 
amplitude that we use is available only for pn energies up to 1.3 GeV, see j3] for details. 

In Fig. [51 we show the four asymmetries as functions of the angle for a fixed missing 
momentum value of p m = 0.4 GeV, and for a fixed P = 35°. The non-relativistic, factor- 
ized PWIA prediction for is oc 1 — 3cos 2 #, which leads to zeros for 9 = 54.7° and 
9 = 125.3°. If we use only the S-wave and D-wave contributions to the ground state wave 
function, and perform the PWIA calculation with the quantization axis along the three- 
momentum transfer, q, we observe exactly this type of angular dependence. The P-wave 
contributions lead to slight deviations from the non-relativistic angular pattern. In the fig- 
ures we show, we have used a quantization axis along the beam, and this rotation obscures the 
original structure of the asymmetry. The angular dependence of A Y ed even in non-relativistic 
PWIA is more complicated than the structure for A d , as the beam- vector asymmetry A y ed is 
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FIG. 5: (Color online) The asymmetries A d (panel (a)), A T d (panel (b)), A\ d (panel (c)) and A^ d 



(panel (d)) for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , p m = 0.4 GeV, and <j) p = 35° are shown 
calculated in PWIA (dotted line), with on-shell FSI (solid line), and with on-shell and off-shell FSI 
(dashed line), as a function of the polar angle of the missing momentum. 



equal to the ratio of helicity-dependent and helicity independent responses. This prevents 
the cancelations of helicity-independent expressions in numerator and denominator that is 
present in the tensor asymmetry A d , and causes its simple angular structure. For A y d and 
A T ed , the non-relativistic result predicts zero for all angles, and this result persists for our 
fully relativistic calculation, for the reasons discussed above. 

Again, we observe the same pattern that was apparent for the momentum distributions: 
for A T d and A^ d , the FSI effects are small for small angles, and become important only 
for larger angles. The differences between on-shell FSI calculations and FSI calculations 
including off-shell FSIs, too, is very small. For A d and A T edl the off-shell FSI effects are more 
pronounced, in particular for A d . Note that when calculating an angular distribution for 
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a fixed missing momentum, we slice through various values of x, and therefore the relative 
importance of the off-shell FSI contributions is different for different angles 9. 

At the four- momentum transfer of 2 GeV 2 , we are limited in the range of polar angles 
9 that we may access, as complete np scattering amplitudes are available only up to 1.3 
GeV. The results for Q 2 = 1 GeV 2 and otherwise identical kinematics are not qualitatively 
different from what we see at smaller angles. 
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FIG. 6: (Color online) The asymmetries A d (panel (a)), A T d (panel (b)), A^ d (panel (c)) and A T ed 
(panel (d)) for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , p m = 0.2 GeV, and 4> p = 35° are shown 
calculated in PWIA (dotted line), with on-shell FSI (solid line), and with on-shell and off-shell FSI 
(dashed line), as a function of the polar angle of the missing momentum. 



We show the asymmetries at a lower missing momentum value, p m = 0.2 GeV, as function 
of the angle in FigEl Overall, it is clear that for the lower missing momentum value, p m = 0.2 
GeV, the influence of FSIs is not that large. As before, the asymmetry A y d shows an off-shell 
FSI result that differs from the on-shell FSI for a larger range of angles, but the effect is 
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much less pronounced than for higher p m . Both A v ed and A T ed show small off-shell FSI effects 
in the region where the asymmetries are large. The off-shell FSI contributions are somewhat 
limited here, as for p m = 0.2 GeV, the maximum kinematically possible x value is 1.3. 

Summarizing, it is interesting to note that the tensor asymmetry A T d and the double 
spin asymmetry A^ d exhibit rather similar behavior, even though they have quite different 
structures: the former depends on the helicity independent terms of the cross section and 
has a tensor (T 20 ) structure, whereas the latter depends on the helicity-dependent terms 
of the cross section, and has a vector structure (T 10 ). Due to invariance under parity and 
time reversal, both responses are non-zero in PWIA, and they show similar structures and 
sensitivity to FSI effects. Their <ft p dependence is similar, too, showing a mirror symmetry 
along P = 180°. In the same way, the target spin asymmetry A d (helicity independent, 
vector) and the tensor-beam asymmetry A T ed (helicity-dependent, tensor) show similar traits: 
they are both zero in PWIA, and are more sensitive even to off-shell FSI effects. Their <j) p 
dependence leads to an inversion of all features above P = 180°. 

3. Contributions from individual parts of the NN scattering amplitude to the FSIs 

In our calculation of the final state interactions, we use the full nucleon-nucleon scatter- 
ing amplitude. There are several ways to decompose and parametrize the NN scattering 
amplitude. It can be parametrized with five terms: a central, spin-independent term, a 
spin-orbit term, and three double-spin flip contributions. It can also be given in terms 
of invariants, using a scalar, vector, tensor, pseudoscalar, and axial term. Some of these 
parametrizations may be more or less useful and enlightening in trying to understand what 
is happening. As we are interested in the effects of target polarization, investigating the 
effects of spin-dependent terms in the FSIs is a logical and interesting step. We separate 
the NN amplitudes into a central term, a single spin-flip (i.e. spin-orbit) term, and three 
double spin-flip terms. 

Fig. [7] shows the contributions of the central, central and single spin-flip, and full FSIs 
to the four asymmetries at Q 2 = 2 GeV 2 and p m = 0.4 GeV as a function of the angle of the 
missing momentum. The tensor asymmetry A T d shows little sensitivity to the details of the 
FSIs, it just shows some minor quantitative changes in the dip and peak region. The beam- 
vector asymmetry, A^ d , is insensitive at lower angles, but shows small changes in magnitude 
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FIG. 7: (Color online) The asymmetries A d (panel (a)), A T d (panel (b)), A\ d (panel (c)) and A^ d 
(panel (d)) for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , p m = 0.4 GeV, and <j) p = 35° are shown 
calculated with on-shell FSI (solid line), without any double spin- flip terms in the on-shell FSI 
(dashed line), and with central on-shell FSI only (dash-double-dotted line), as a function of the 
angle 9 m of the missing momentum. 

at larger angles. In both cases, there is no shape change when the different parts of the FSIs 
are added. This changes when considering the target-spin asymmetry A d and the tensor- 
beam asymmetry A T ed . For these asymmetries, the shape is quite different when only the 
central part of the FSIs is included. The central FSI result for Aj d is rather small and even 
takes some negative values in a shallow dip around 9 m 45°. Once the single spin-flip FSI is 
included, the asymmetry changes and shows a steep rise with a shallow shoulder at larger 
angles. With the inclusion of the double-spin flip FSIs, the magnitude of the asymmetry 
increases a bit at larger angles. 

For A d , the influence of the spin- dependent FSIs is most pronounced: while the asym- 



19 



metry is very small and changes sign twice with central FSIs, the inclusion of the single 
spin-flip term leads to an asymmetry that is similar in shape, albeit a bit larger than with 
central FSIs, and of opposite sign. The double spin-flip terms completely change the shape 
of the asymmetry, leading to a pronounced peak and a much larger maximum value. Here, 
for A d , the effect of the double spin-flip FSIs is most pronounced and most relevant. 




FIG. 8: (Color online) The asymmetries A d (panel (a)), A d (panel (b)), A^ d (panel (c)) and 
A T ed (panel (d)) for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , x = 1, and 4> p = 35° are shown 
calculated with the full on-shell FSI (solid line), without any double spin- flip terms in the on-shell 
FSI (dashed line), and with central on-shell FSI only (dash-double-dotted line) as a function of the 
missing momentum. 

Fig. [8] shows the contributions of the central, central and single spin-flip, and full FSIs 
to the four asymmetries at Q 2 = 2GeV 2 and function of the missing momentum. 

While the angular distributions for p m = 0.4 GeV shown in Fig. [7] do not show a very 
pronounced effect of the double spin-flip terms on A d and A^ d , the momentum distribution 
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for A T d shows that especially for larger missing momenta, the spin-dependent FSIs are very 
relevant. Starting for p m = 0.4 GeV, the results without the full spin-dependence deviate 
significantly from the full FSI result, and for p m > 0.5 GeV, an interesting inversion happens: 
the result for central FSI only is closer to - but still far off - the full FSI result than the 
calculation without double spin-flip. This indicates that interference effects are relevant for 
A T d in this kinematic region. For A d , a similar picture emerges for larger missing momentum: 
for p m > 0.4 GeV, central FSI only results are above the full result, while the no double spin- 
flip result is below it. For this asymmetry, for p m > 0.1 GeV, all types of spin- dependent 
FSI are very important. The tensor-beam asymmetry, A^ d , shows that while the double 
spin-flip terms have only a small effect, the single spin-flip term, i.e. the spin-orbit term, 
gives a huge contribution. 

Overall, the tensor asymmetry A T d and the double spin asymmetry A^ d exhibit rather 
similar behavior, showing some quantitative and no large qualitative dependence on spin- 
dependent FSIs. The target-spin asymmetry A d and the tensor-beam asymmetry A T ed show 
large qualitative and quantitative sensitivity to spin-dependent FSIs, each in a different way. 

In a previous paper 7|] dealing with unpolarized observables, we investigated the influence 
of the different invariant amplitudes of the NN amplitude parametrization by calculating 
the FSIs with only one of the invariant amplitudes. For the unpolarized case, we found that 
the role of interference is huge, and that there is no single dominant amplitude. For the 
asymmetries, we find that for small angles, the pseudoscalar amplitude seems to be very 
close to all asymmetries except for Aj d , but this behavior is confined to 9 < 30°. Deviations 
beyond that are significant, in particular for A d . The results show that overall, there are 
many relevant interference effects, and no single part of the NN amplitude is dominant. 



4- Influence of the D wave 



A question often discussed is the influence of correlations in the nuclear ground state, 
and in the case of the deuteron, the role played by D-wave - and P-wave - admixtures. Due 
to the rather different predictions of various non-relativistic iViV models for the D-wave 
content, the hunt for observables sensitive to this part of the wave function has been going 
on for a long time. It should be noted that on theoretical grounds the attempts to extract 
the D-wave contribution to the deuteron bound state is ill considered. The wave function 
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is not an observable and unitary transformations can change the D-state contribution while 
leaving the matrix elements unchanged. Thus an actual observable contains information 
about initial and final states, as well as on the current operator, with the various quantities 
changing with unitary transformation and one cannot be uniquely separated from the others. 

In our calculation, due to the normalization of the ground state wave function, there 
are some issues with directly isolating the D-wave contribution. Just in order to give an 
impression of the influence of the D-wave contribution on the asymmetries we study here, we 
have simply switched off the D-wave contributions, without changing the normalizations. In 
our relativistic calculation, there is also a P-wave contribution present. We study its effect, 
too. In general, P-wave effects are expected to be very small. 




FIG. 9: (Color online) The asymmetry Aj, for a beam energy of 5.5 GeV, Q 2 = 2 GeV 2 , and 
p m = 0.4 GeV is shown calculated (a) with on-shell FSI, (b) without a D-wave contribution to the 
ground state, (c) without a P-wave contribution to the ground state, and (d) with only an S-wave 
contribution to the ground state, as a function of polar angle of the missing momentum 9 m and 
the proton's azimuthal angle. 

While the other asymmetries also show significant dependence on the D-wave, we will 
focus here for brevity on the effects of the D-wave and P-wave contributions on the tensor 
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asymmetry A d . As expected, the difference between the calculation with the full ground- 
state wave function (a) and the calculation without the D-wave (b) is large: a prominent dip 
is turned into a peak, and the maximum values reached change. In non-relativistic PWIA, 
this asymmetry would be zero without the D-wave, but with FSI - even just central FSI - 
the tensor asymmetry acquires a nonzero value, as the relative position of the neutron and 
the knocked-out proton matter for the strength of the FSI that is experienced. 

Performing a calculation without the P-wave contribution (c) does not lead to any signif- 
icant changes, the peak heights vary a little, but there are no qualitative changes. Panel (d) 
shows the results for just the S-wave part of the wave function. Here, the missing P-wave 
contribution - still present in the top right panel without the D-wave - leads to a somewhat 
different shape and an increased magnitude for the dip structure at lower <p p values. It 
is interesting to note that in PWIA, if we switch off the D-wave contribution, the tensor 
asymmetry is small but still non-zero due to the P-wave contributions. 

IV. SUMMARY AND OUTLOOK 

In this paper, we have presented a formalism for the calculation of responses and asym- 
metries for polarized deuteron targets. We have shown how to evaluate these observables 
in different reference frames, and for different polarization axes. Symmetries of the current 
matrix elements were pointed out, and together with the behavior under parity and time 
reversal transformations, exploited to show that two of the asymmetries we discuss, the 
target-spin asymmetry A v d and the tensor-beam asymmetry Aj d , vanish in PWIA. 

We performed a relativistic calculation of various asymmetries accessible with a polar- 
ized deuteron target. We have included a full FSI calculation, with on-shell and off-shell 
contributions, using experimental data on the pn scattering amplitude as input. Final state 
interactions are very relevant for all asymmetries in most kinematics. Two of the asymme- 
tries vanish in PWIA, and are therefore more sensitive to FSI effects, even to the off-shell 
FSI contributions. An important result of our paper is that even in the region of the quasi- 
elastic ridge, x — 1, the influence of FSIs on the asymmetries is large, and a straightforward 
extraction of D-wave properties from measured data will not be possible. This is true even 
though the influence of the D-wave on the asymmetries is large, as commonly assumed. The 
influence of the P-waves, a purely relativistic phenomenon, is generally small, unless we 
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consider situations where the D-wave is switched off. 

One interesting and conspicuous feature of the asymmetries is the fact that the target- 
spin asymmetry A v d and the tensor-beam asymmetry A T ed have very similar properties - they 
vanish in PWIA, have similar sensitivity to FSIs, and a similar dependence on <p p - just as 
the tensor asymmetry A T d and the beam-vector asymmetry A^ d show similar properties in 
these respects. We have shown that this can be understood in terms of their behavior under 
parity and time reversal in PWIA. 

We have tested the sensitivity of our results to the different parts of the FSIs. As expected, 
spin-dependent FSIs are relevant, and depending on the kinematics and observable, even the 
double spin-flip terms are extremely important. 

Our calculation has been performed in impulse approximation, i.e. assuming that the 
detected proton is the nucleon that interacted with the photon initially. Contributions from 
the Born term, where the photon interacts with the neutron, will in general be small for 
most observables in most kinematics, but they may contribute noticeably for larger missing 
momenta (p m > 0.6 GeV). 

Next, we plan to perform calculations for an unpolarized deuteron target and a polarized 
ejected nucleon. 
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APPENDIX A: THE DENSITY MATRIX 

Consider an object with total angular momentum j and projection m. An arbitrary state 
of angular momentum j can be written as 




(Al) 



m=-j 



where normalization of the state requires that 



E H 



1 . 



(A2) 



m=-j 
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The expectation value of some operator A for this state is given by 

3 3 j j 

< A > i= c in* A \jm) 4 = E E & m 'l A l jm > 44* • (A3) 

m'=j rn=—j m'=j m=—j 

Any attempt to polarize a target consisting of a collection of these objects by applying 
magnetic fields will in general not produce a single state such as that described above, but 
will consist of a statistical ensemble of such states with probabilities Pi such that 

$> = !■ (A4) 

i 

The statistical average of the expectation value of operator A is then given by 

< a >= J2 Pi < a >i= z E & m 'i 1 i-?' m > E « c - ; • ( A5 ) 

i m'=—j m=—j i 

Defining the density matrix 

Pmm' = Z] C mPi C ln' (A6) 
i 

and 

A m / m = (j,m'\A\jm) (A7) 

the average expectation value of A can be written as 

3 j 

< i >= E E A m'mPmm' = Tr(Ap) (A8) 
m'=—j m=—j 

where A and p are matrix representation of A and the density matrix in the subspace of 
total angular momentum j. 
From flA6l) . 

Pmm' = ^ y C m'Pi C m = Pm'm (A9) 
i 

or in matrix form 

pt = p . (A10) 

So the density matrix is hermitian. Also, 

t Kp) = E E«4* = E E I4 2p = E p * = !■ ( An ) 

m=—j i i m=—j i 
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A further constraint on density matrix is given by 

Tr(p 2 )<(Tr(p)) 2 = l. (A12) 

It is often convenient to express the density matrix for angular momentum j in terms of 
spherical tensor operators such that 

1 23 J 

p = 27TI £ £ T *^ fjM (A13) 

J ' J=o M=-.J 

where tjm is an irreducible spherical tensor operator of rank J and projection M and the 
T JM are complex coefficients that describe the average polarization of the target. These are 
defined such that 

T* JM = (-1) M Tj-m (A14) 

and 

t'jm = (-i) M rj- M ■ (A15) 

Since we are concerned with a polarized deuteron target in this paper we will now confine 
the argument to the case j = 1. In this case matrix elements of the density operator are 
given by 

1 2 J 

Pav = (1A|/)|1A') = -^ T* jm {1\\t jm \1\') . (A16) 

J=0 M=-J 



This can be written in matrix form as 

2 

3 



P=lt, E T JmTjm (A17) 

J=0 M=-J 

If we choose normalizations such that 

T o = 1 (A18) 

and 

(1 ||fj|| 1) = V3V2J+1, (A19) 
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the matrices tjm are 



no 



no 
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72 



^1 ^ 

1 
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0-1 

^10 0^ 

0-2 
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Til 



ni 



0-10 
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f -1 o\ 
1 




T~22 



and the remaining matrices can be obtained from 



' JM 



:-i) M n_M 



These matrices have the properties 



Tr(r JM ) = 



and 



So, 



Tr ( T j'M' T JM) = 3 5j/j5 



M'M 



Tr(rt M p) = Tj M . 
The constraint given in (IA12I) requires that 



2 J 



< 1 



J=l M=-J 

Using the matrices defined by ( 1A20I) . we can write 



j=i 



TjqTjq + (T} M T JM + Tj_ m Tj- M ) 



M=l 
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The last term of this can be rewritten using (1A14I) to give 



^ ( t jm t jm + Tj-mT j-m) = 7^ (Tjm t jm + (-1) M T jm tj-m) 



M=l 



M=l 
J 



[^{Tjm) (tjm + {-1) M Tj- M ) 



M=l 



+3(T JM )(-z) (t jm -{-1) m Tj„ m )] 

J 

(®(Tj M )rf M + %(T jm )tJ m ) (A27) 



M=l 



where 



r 5R 

JM 



t 

JM 



(A28) 



and 



' JM 



{tjm - (-1) m tj_ m ) = -i (tjm - t\ m 



The orthogonality relations for the tjm can be used to show that 
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